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SYNOPSIS 


Five  additional  embankment  models  have  been  tested  to  complete  a 
research  project  on  the  verification  of  the  Nawmark  (1965)  sliding  block 
analysis  to  predict  permanent  displacements  of  earth  dams  subjected  to 
earthquake  loading.  The  tests  were  conducted  on  a  new  "Bumpy  Road" 
shaking  table  aboard  the  Cambridge  University  Geotechnical  Centrifuge. 

Clay  embankment  models  of  triangular  cross  section  with  1:2  slopes  that 
represent  prototypes  of  6.9  m  height  were  each  submitted  to  between  three 
and  six  earthquakes.  The  sliding  block  calculations  consistently  under¬ 
estimated  displacements;  however,  relative  errors  in  displacements  were 
smaller  for  large  than  small  earthquakes .  The  yield  of  clay  at  stresses 
below  the  shear  strength,  dynamic  magnification  of  base  acceleration  and 
strain  softening  during  shaking  are  cited  as  factors  contributing  to 
unconservative  predictions  of  displacement. 

Another,  unexpected  mode  of  embankment  failure  was  observed:  strain 
softening  and  pore  pressure  generation  during  one  test  caused  a  slope 
failure  3.5  s  (two  hours  in  prototype  terms)  after  the  model  earthquake. 

The  displacement  during  this  slung;  was  about  2.8  m  in  prototype  terms  - 
on  the  order  of  ten  times  the  displacements  which  occurred  during  any  actual 


earthquake 
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2.  INTRODUCTION 

The  purpose  of  this  research  was  to  test  the  validity  of  using 
Newmark  sliding  bl jck  calculations  to  predict  permanent  deformations  of  clay 
embankments  during  earthquakes.  In  view  of  the  unpredictability, 
unrepeatability  and  Infrequency  of  real,  large  earthquakes ,  It  Is  valuable 
to  conduct  model  tests  to  improve  our  understanding  of  the  behaviour  of 
earth  structures  during  base  shaking. 

Soil  models  at  one  gravity  can  help  to  supplement  field  data  but 
suffer  from  the  drawbacks  of  using  either  very  large  (difficult  to  prepare) 
models  or  small  models  at  very  low  stress  levels.  Centrifuge  modelling 
allows  the  researcher  to  produce  small  (easy  to  prepare)  models  at  the 
correct  stress  levels.  Hence,  it  is  a  useful  tool  for  producing  experimental 
data  to  test,  verify  and  develop  our  understanding  of  the  dynamic  behaviour 
of  soil. 

It  was  proposed  to  conduct  a  series  of  model  tests  utilizing  an 
existing  shaking  table  aboard  the  Cambridge  Centrifuge.  It  was  proposed 
that  by  varying  the  clay  shear  strength,  the  centrifuge  speed,  and  the 
earthquake  magnitude,  the  deformation  during  a  range  of  combinations  of 
yield  acceleration  and  base  acceleration  could  be  measured  for  comparison 
with  sliding  block  theory. 

For  the  five  initial  tests  reported  in  the  first  progress  report 
(M.  Phil,  thesis,  B.Ii.  Kutter)  it  was  attempted  to  keep  approximately  the 
same  shear  strength  in  all  models  and  only  to  vary  the  centrifuge  acceleration 
and  the  earthquake  acceleration.  These  tests  gave  some  good  results  but 
also  revealed  problems  and  suggested  possible  improvements  in  the  test 
procedure  as  well  as  making  clear  the  difficulty  of  carrying  out  some 
tests  suggested  in  the  proposal: 

(1)  It  is  impractical  to  test  models  with  yield  accelerations  of  5%  of 
the  centrifuge  accelerations  or  less.  This  would  imply  static  factors  of 


safety  on  the  order  of  1.1  which  besides  being  atypical  of  real  dams,  could 
result  in  ruining  a  model  due  to  creep,  small  errors  in  centrifuge 
speed  or  clay  shear  strength. 

(2)  The  geometry  of  models  used  for  the  initial  tests  was  undesirable  in 
two  respects.  Firstly,  the  embankment  slopes  were  steeper  than  those 
typical  of  earth  dams.  Secondly,  there  was  only  one  sloping  face,  and  the 
clay  was  confined  on  the  other  side  by  the  wall  of  the  shaking  table. 

Although  this  wall  was  well  clear  of  "static"  (non-earthquake)  failures, 
earthquake  shear  zones  were  found  to  be  more  reclined  and  did 
actually  intersect  this  clay  wall  boundary. 

(3)  Using  the  spring  mass  resonant  shaking  table  on  clay  specimens  was 
difficult  since  only  one  earthquake  could  be  released  per  centrifuge  flight. 
In  addition,  the  catch  mechanism  could  not  be  cocked  until  just  before  a 
model  earthquake  since  premature  firing  of  the  earthquake  would  occur  in  a 
matter  of  tens  of  minutes  after  cocking.  These  conditions  made  it  impossible 
to  achieve  complete  pore  pressure  equilibrium  before  the  earthquakes  and 
also  necessitated  numerous  stop-start  centrifuge  cycles  in  the  course  of 
model  testing. 

For  these  and  other  reasons,  a  new  actuating  system  was  developed  for 
additional  tests.  The  new  apparatus  afforded  much  more  versatility  but 
also  introduced  one  different  test  restriction.  It  was  not  capable  of 
imposing  earthquakes  greater  than  40%  of  the  centrifuge  acceleration.  This 
precluded  the  possibility  of  performing  the  large  earthquakes  originally 
proposed  at  the  low  centrifuge  accelerations. 

The  new  apparatus  allowed  construction  of  larger  embankment  models 
with  two  sloping  faces,  each  with  1:2  slopes.  In  order  to  limit  the 
number  of  new  variables  it  was  considered  desirable  in  the  additional  tests 
to  again  hold  the  clay  shear  strength  constant.  Since  the  magnitude  of 
the  earthquakes,  the  centrifuge  acceleration,  and  the  frequency  of  vibration 
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were  the  easiest  and  most  accurately  adjustable  variables,  these  were  the 
only  parameters  that  were  varied  in  the  additional  tests. 
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THE  BUMPY  ROAD  APPARATUS 


The  apparatus  was  conceived  by  P.W.  Turner  and  the  design  was  carried 
out  jointly  by  P.W.  Turner  and  B.L.  Kutter.  It  operates  by  forcing  a 
wheel  on  the  centrifuge  arm  to  follow  a  curved  steel  track  that  is  attached 
to  the  wall  of  the  centrifuge  pit.  The  wheel  vibrations  are  then  transmitted 
through  a  crankshaft  and  slider  mechanism  to  vibrate  the  soil  model 
tangentially  relative  to  the  centrifuge.  Several  earthquake  tests  on 
board  the  Cambridge  centrifuge  have  been  conducted  using  a  spring  mass 
resonant  actuator  designed  by  Morris  (1979) ,  however  this  apparatus  was 
limited  in  sample  dimension,  and  only  a  single  earthquake  could  be  released 
during  a  centrifuge  flight. 

The  Bumpy  Road  actuator  provides  improved  dynamic  motions  to 
centrifuge  samples  with  the  following  features: 

*  20  g  "earthquake"  accelerations. 

*  Vibration  frequencies  up  to  about  150  Hz. 

*  Capacity  for  135  kg  models. 

*  Centrifuge  accelerations  up  to  100  g. 

*  A  soil  container  with  approximate  dimensions  of  0.9  x  0.48  x  0.22  m 
(length:width:height)  to  accommodate  embankment  dam  models . 

*  Provision  for  viewing  a  cross  section  of  the  models  during 
flight  by  using  flash  photography. 

*  Ability  to  release  several  earthquakes  during  flight  without 
stopping  the  centrifuge. 

*  Capability  of  in-flight  variability  of  earthquake  magnitude. 

*  Provision  for  changing  earthquake  waveform  by  changing  the  track 
on  the  centrifuge  wall. 

Two  schematic  views  of  the  system  are  shown  in  Figure  3.1.  The 
track  (1)  extends  over  about  1/3  of  the  circumference  of  the  centrifuge  pit 
and  is  machined  to  a  prescribed  waveform  by  a  numerically  controlled  milling 
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machine.  As  the  wheel  (2)  follows  the  track  the  crank  (3)  and  shaft  (4) 
are  rotated  dynamically.  The  wheel  Is  held  against  the  track  by  the  double 
acting  pneumatic  jack  (5).  The  vibration  take  off  slider  (6)  can  be  moved  to 
an  eccentric  position  relative  to  the  shaft  by  an  electric  motor  and  screw 
jack  (7)  to  provide  any  desired  ratio  of  package  displacement  to  wheel 
displacement  between  0.0  and  0.2.  The  slider  itself  imparts  the  tangential 
dynamic  loads  to  am  actuating  beam  (8)  which  rests  on  linear  roller 
bearings  (9) .  The  beam  is  bolted  to  a  rack  (10)  which  engages  with  a 
matching  rack  on  the  swinging  platform  (11) .  In  plan  view  looking  vertically 
into  the  centrifuge  pit,  the  model  dam  (12)  is  enclosed  on  the  ends  by 
aluminium  alloy  walls  (13).  Slumping  soil  from  a  model  can  be  allowed  to  fall 
into  the  overflow  chambers  (14)  to  prevent  the  build  up  of  soil  against 
the  end  walls.  The  model  is  viewed  in  cross  section  through  polycarbonate 
windows. 

The  platform  itself  is  suspended  by  tangentially  flexible  straps  (16) 
from  torsion  bars  (17) .  As  the  centrifuge  acceleration  is  increased  from 
zero  the  package  swings  up  until  the  matching  racks  become  aligned.  The 
torsion  bus  allow  the  package  to  move  radially  outward  relative  to  the 
centrifuge  arm  (18)  until  the  matching  racks  engage  causing  the  additional 
centrifugal  load  to  be  transferred  to  the  end  of  the  arm  through  the 
linear  roller  bearings. 

The  tracks  are  fixed  to  the  wall  by  fifteen  brackets  (19)  which  are 
in  turn  fixed  to  the  wall  (20)  by  Rawl  Bolts.  The  bracket  spacing  is  ten 
degrees  of  arc  of  the  centrifuge. 

The  pneumatic  system,  extremely  simplified  in  figure  3.1  is 
represented  by  two  reservoirs  (21)  and  two  solenoid  valves  (22)  and  (23) . 

Part  of  the  pneumatic  system  is  actually  mounted  at  the  centrifuge  axis . 

There  is  an  electronic  triggering  system  which  synchronises  the  opening 
and  closing  of  the  solenoid  valves  with  the  position  of  the  arm  of  the 


centrifuge  relative  to  the  track.  So,  at  the  push  of  a  button  in  the 
console  room,  valve  (23)  is  closed  and  (22)  is  opened,  pulling  the  wheel 
out  to  follow  the  track  for  a  single  revolution.  Then  (23)  is  opened  and 
(22)  is  closed  so  that  the  wheel  is  pushed  in  to  miss  the  track  on 
subsequent  passes.  Photographs  of  the  assembled  apparatus  from  various 
points  are  shown  in  Figures  3.2a  to  3. 2d. 
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INITIAL  TESTS:  SUMMARY  OF  RESULTS 
A  relationship  between  normalised  displacement  and  the  ratio  of 
yield  accelerations  to  base  accelerations  km  was  developed  for  sinusoidal 
base  acceleration  and  is  shown  in  figure  4.1.  The  predictions  of  displacements 
using  this  relationship  for  the  initial  tests  on  the  spring  mass  actuator 
consistently  underestimated  measured  displacements.  However,  at  the 
conclusion  of  the  initial  tests  it  was  unclear  whether  the  errors  in 
prediction  were  a  result  of  deficiencies  in  the  sliding  block  calculations 
or  deficiencies  in  the  soil  models.  Evidence  from  radiographs  indicated 
that  the  assumption  that  a  distinct  rupture  plane  forms  during  sliding 
was  not  accurate.  But  the  importance  of  this  inaccuracy  could  not  be 
evaluated. 

The  radiograph  for  model  BK-5  shown  in  figure  4.2  illustrates  the 
difference  between  dynamic  deformations  (which  occur  in  a  wide  band)  and 
static  failure  (which  occurs  on  a  rupture  surface) . 

Sample  BK-5  was  first  subjected  to  earthquake  deformation  which 
occurred  on  the  wide  band  that  intersects  the  left  side  of  the  photograph. 

After  the  earthquake  the  centrifuge  acceleration  was  increased  until  a 
"static"  slip  occurred  on  the  very  distinct  circular  arc.  that  intersected 
the  top  surface  of  the  specimen. 

There  are  many  soil  properties  which  influence  the  tendency  to  form 
rupture  surfaces .  But,  in  this  case  the  reason  for  the  differences  between 
dynamic  and  static  mechanisms  is:  for  static  (constant)  loading,  failure 
occurs  when  the  factor  of  safety  is  one,  so  that  soil  only  on  the  most 
critical  slip  surface  is  subjected  to  failure  stresses  and  deformation  is 
confined  to  this  surface;  in  contrast,  dynamic  loading  can  produce  factors 
of  safety  much  less  than  one  for  a  short  time.  During  this  time,  deformation 
occurs  on  a  band  of  slip  surfaces  which  all  have  factors  of  safety  less 


than  one. 


In  spite  of  the  discrepancy  between  assumed  and  observed  mechanisms , 
sliding  block  displacement  predictions  seemed  reasonable.  However,  the 
fact  that  predictions  cure  highly  sensitive  to  the  assumed  value  of  soil 
strength  had  to  be  taken  into  account.  To  demonstrate  this  point,  consider 
an  example.  A  sliding  block  with  a  yield  acceleration  of  k^  g  of  0.1  g  rests 
on  an  inclinedl:2  cohesive  slope  (see  figure  4.3).  A  single  horizontal 
sinusoidal  acceleration  pulse  of  amplitude  k^  g  at  a  frequency  of  1  Hz  is 
applied  to  the  base  of  the  model.  Now  consider  the  error  in  displacement 
predictions  during  the  single  acceleration  pulse  if  the  predictions  are  based 
on  a  5%  overestimate  of  shear  strength  (errors  of  5%  in  the  determination 
of  shear  strength  would  not  be  unusual) .  For  a  large  sine  pulse  with 
ky/k^  *0.2  one  would  predict  a  displacement  of  400  mm  and  would  observe 
460  am,  a  15%  unconservative  error.  For  one  small  sine  pulse  with  k^/k^  *0.9 
a  displacement  of  1.0  mm  would  be  predicted  and  17  mm  would  be  observed, 
a  1600%  unconservative  error.  Relative  error  tends  to  be  larger  for  small 
pulses  that  barely  exceed  the  yield  acceleration.  This  fact  makes  it 
difficult  to  test  and  compare  analytical  techniques  based  on  existing  field 
evidence  of  measured  displacements  during  relatively  small  earthquakes . 

The  sensitivity  of  dynamic  displacements  to  strength  is  the  nature  of  the 
problem,  not  of  the  analytical  technique  and  therefore  should  not  necessarily 
be  considered  a  drawback  of  sliding  block  type  calculations. 

In  conclusion,  the  relative  importance  of  all  assumptions  of  an 
analytical  technique  should  be  kept  in  mind.  The  assumption  that  a  distinct 
failure  plane  forms,  albeit  inaccurate,  may  be  less  significant  them 
inaccurate  characterisation  of  soil  strength.  Centrifuge  modelling  allows 
one  to  conduct  a  repeatable  series  of  tests  helping  us  to  keep  these 
considerations  in  the  proper  perspective. 


5. 


BUMPY  ROAD  TESTS:  PREPARATION  AND  PROCEDURE 


All  of  the  nodal  earthquake  teats  reported  here  were  conducted  on 
anbanknents  of  triangular  cross  section,  with  1:2  slopas.  The  models  were 
tr lamed  from  blocks  of  kaolin  clay  consolidated  in  a  rectangular  tub  to  a 
vertical  stress  of  124  kPa.  The  blocks  of  clay  were  extruded  (see  figure 
5.1)  and  then  trimmed  to  the  desired  profile  (see  figure  5.2).  Sample 
dimensions  and  typical  transducer  locations  are  shown  in  figure  5.3  and 
5.4.  Some  pore  water  pressure  transducers  were  suspended  in  the  clay  at 
the  time  of  pouring  the  slurry  into  the  conso Udometer .  Other  pore  water 
pressure  transducers  were  inserted  along  with  accelerometers  and  LVDT 
resting  pads  after  trimming  the  sample. 

Lead  powder  threads  were  inserted  into  the  specimens  before  testing 
to  allow  post  test  x-radiography  to  reveal  the  interned,  deformation  mechanisms 
caused  by  earthquakes.  Post  test  radiographs  are  shown  in  figures  5.5a  to 
5.5e.  Silver  coated  nylon  markers  were  placed  in  grids  at  10  mm  spacing 
to  enable  in  flight  photographs  to  reveal  the  deformation  patterns  at  the 
specimen  boundary.  Figures  5.6a  to  5.8b  are  examples  of  photographs  which 
were  taken  at  various  stages  during  the  centrifuge  tests.  For  more 
information  about  specimen  preparation  techniques  refer  to  the  first  progress 
report  (M.  Phil,  thesis,  B.L.  Kutter)  written  for  this  proposal  or 
B.L.  Kutter 1 s  forthcoming  Ph.D.  thesis. 

Post  test  photographs  of  model  C  are  shown  in  figures  5.9a  to  5.9b, 
not  only  because  of  the  dramatic  failure,  but  to  give  the  reader  a  feel  for 
the  appearance  of  an  instrumented  model. 

The  basic  test  procedure  was ■ to  bring  the  centrifuge  up  to  the  test 
speed  then  fill  the  reservoirs  to  the  desired  height  and  allow  the  pore 
water  pressures  to  come  into  equilibrium.  After  each  earthquake  was 
released  the  sample  was  allowed  to  consolidate  before  the  next  earthquake. 

The  magnitude  of  successive  earthquakes  was  adjusted  if  desired.  It  is 
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Important  to  allow  the  water  pressure  in  a  model  to  come  into  equilibrium 
if  accurate  estimates  of  shear  strength  are  to  be  deduced. 

Table  I  summarises  the  characteristics  of  the  5  models  tested  on  the 
Bumpy  Road.  Most  of  the  model  tests  simulated  a  6.9  m  high  embankment 
shaken  at  a  prototype  frequency  of  1.52  Hz  (tests  A,  B,  C,  D) .  Models  A 
and  B  were  tested  at  80  g  and  were  therefore  2.25  times  smaller  than 
Models  C  and  D  which  were  tested  at  35.6  g.  Because  of  the  scaling 
relationship  for  dynamic  centrifuge  modelling  (see  first  progress  report) 
models  A  and  B  had  to  be  shaken  at  a  higher  frequency  than  C  and  D,  and 
were  therefore  conducted  with  a'  bumpy  road  track  that  had  a  smaller  wave¬ 
length.  Model  E  was  a  slightly  larger  model  than  C  end  D  and  was  tested  at 
35.6  g  on  the  same  wavelength  track  to  test  the  effect  of  varying  the 
prototype  frequency  from  1.52  to  2.28  Hz. 

Another  variation  from  the  standard  frequency  and  model  scale  were 
tests  AXV,  the  fourth  earthquake  on  Model  A  and  test  CV,  the  fifth  earth¬ 
quake  on  Model  C.  These  tests  were  conducted  after  rapidly  increasing 
the  centrifuge  to  1.25  times  the  consolidation  acceleration  and  releasing 
an  earthquake  before  the  sample  had  a  chance  to  consolidate. 


6. 


BUMPY  ROAD  TESTS:  SOME  RESULTS  AND  ANALYSIS 


6.1  INTRODUCTION 

This  section  begins  with  a  description  of  data  collected  during  two 
model  earthquakes.  The  first  one  to  be  discussed  is  test  BI,  the  first 
earthquake  on  one  of  the  small  models  at  80  g.  The  second  to  be  discussed 
is  from  test  DX,  the  first  earthquake  on  one  of  the  large  models  tested 
at  35.6  g. 

The  section  continues  with  a  discussion  of  the  determination  of  the 
soil  shear  strength  and  the  yield  acceleration.  Using  these  parameters, 
sliding  block  displacements  were  calculated  and  compared  with  measured 
displacements.  Reasons  for  the  underprediction  of  displacements  by  the 
theory  are  discussed. 

The  last  put  of  this  section  discusses  the  delayed  failure  of 
model  C,  approximately  3.5  seconds  after  earthquake  CV. 

Data  picked  off  magnetic  tape  records  or  digital  data  logger  records 
is  summarised  in  Table  II.  Base  and  crest  peak  filtered  acceleration, 
net  crest  settlement  and  lateral  displacement  of  the  slope  faces  at  mid- 
height  are  listed  for  every  earthquake. 

Before  beginning  this  analysis  it  is  useful  to  discuss  the  repeat¬ 
ability  of  the  model  tests,  and  to  demonstrate  the  agreement  between  model 
tests  at  different  scales.  Models  A  and  B  were  small  models  and  were  shaken 
at  80  g  centrifuge  acceleration.  Models  C  and  D  were  2.25  times  larger 
and  were  tested  at  35.6  g.  The  bumpy  road  track  was  changed  after  model 
tests  A  and  B  so  that  the  frequency  of  vibrations  on  all- four  models  was 
scaled  according  to  the  dynamic  centrifuge  modelling  laws.  All  four  models 
represented  6.9  m  high  embankments  that  were  shaken  at  1.52  Hz.  The 
response  of  all  four  models  should  then  be  similar  -  since  they  represent 
the  same  prototype,  all  aspects  of  the  dynamic  behaviour  should  be  related 
by  the  scaling  laws.  Testing  the  scaling  laws  by  constructing  models  of 
different  scales  that  represent  the  same  prototype  is  called  "modelling  of 
models". 


One  aspect  of  the  agreement  between  the  models  of  different  scale 
has  been  Investigated  and  is  illustrated  in  figure  6.1.  The  prototype  crest 
settlement  is  plotted  as  a  function  of  peak  base  acceleration.  It  is  seen 
that  all  four  models  that  were  shaken  at  a  prototype  frequency  of  1.52  Hz 
and  had  the  same  prototype  height  fall  on  the  same  curve.  Despite  the 
fact  that  the  earthquake  amplitudes  were  not  uniform,  and  were  noticeably 
different  for  the  tests  at  35.6  g  and  80  g,  the  agreement  is  very 
encouraging. 

6.2  DESCRIPTION.  OF  TYPICAL  TEST  DATA 

Examples  of  the  data  collected  during  the  model  earthquakes  are 
shown  in  figures  6.2  and  6.3. 

Figure  6.2  is  the  result  of  the  first  earthquake  on  model  B  (see 
figure  5.3  for  transducer  locations).  The  three  base  accelerations  recorded 
in  figure  6.2a  all  show  very  similar  wave  forms  and  amplitudes .  The 
acceleration  at  mid-height  in  the  clay  shows  some  evidence  of  high  frequency 
vibration  during  the  initial  smaller  amplitude  pulses  and  shows  an  attenuation 
of  the  sharp  peaks  of  acceleration  during  the  larger  amplitude  pulses.  The 
horizontal  crest  acceleration  is  of  a  greater  amplitude  them  that  for  the 
acceleration  at  mid-height,  it  also  exhibits  a  more  pronounced  high 
frequency  component  during  the  initial  earthquake  pulses.  This  high 
frequency  component  is  due  to  resonance  of  the  embankment.  During  the 
large  pulses,  the  crest  acceleration  is  rounded  or  chopped  off  at  the  peaks 
due  to  yielding  of  the  embankment.  This  yielding  causes  a  reduction  of  the 
effective  soil  modulus  and  an  increase  of  damping,  consequently,  the  high 
frequency  ringing  of  the  embankment  disappears.  The  vertical  acceleration 
at  the  crest  of  the  embankment  also  exhibits  a  high  frequency  response 
during  the  initial  pulses  of  the  earthquake,  but  again  this  high  frequency 
response  is  lost  during  the  large  acceleration  pulses  while  the  embankment 
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yields  and  acceleration  pulses  of  a  frequency  similar  to  the  basic  driving 
frequency  appear. 

Figure  6.2b  shows  the  crest  settlement  and  the  horizontal  displacement 
of  the  two  slope  faces  measured  half  way  up.  the  slope.  The  pore  pressure 
transduces  at  (11)  and  (12)  H  are  not  responding  to  the  dynamic  pressures , 
probably  due  to  air  coming  out  of  solution  from  the  water  while  the 
specimen  sat  under  40  kPa  pore  water  suctions  during  preparation.  However, 
the  transducer  at  (13)  height  is  responding  to  the  dynamic  loading  and  displays 
a  generation  of  positive  pore  pressure. 

Figure  6.3  shows  the  acceleration  and  displacement  response  of  the 
first  earthquake  on  model  D  at  35.6  g  (see  figure  5.4  for  transducer 
locations) .  The  amplitude  of  the  base  acceleration  is  more  uniform  than 
that  for  the  tests  at  higher  frequency  (tests  A  and  B) .  As  for  the  earth¬ 
quakes  on  the  small  models,  we  see  some  higher  frequency  response  of  the 
accelerometers  in  the  embankment  during  the  first  pulse.  The  crest 
accelerations  have  flattened  peaks  due  to  yielding  of  the  clay. 

All  of  the  displacement  transducers  reveal  an  acceleration  in  the 
rate  of  deformation  at  about  the  middle  of  the  earthquake  despite  the  fact 
that  the  base  accelerations  are  actually  reducing  in  amplitude  throughout 
the  earthquake.  This  is  pechaps  caused  by  softening  of  the  clay  during  the 
quake' due  to  pore  water  pressure  generation  and/or  remoulding.  Figure  6.3 
shows  the  pore  pressure  response  at  three  locations  in  the  embankment.  The 
transducer  (10)  (on  figure  5.4)  near  the  bottom  of  the  embankment  shows  a 
positive  pore  pressure  generation  of  about  25  kPa  compared  to  a  total  over¬ 
burden  of  about  84  kPa.  The  transducer  at  (10)  also  exhibits  a  frequency 
response  of  twice  the  bate  shaking  frequency  -  possibly  due  to  the  soil 
going  through  a  contraction/dilation  cycle  every  time  the  shear  stress  is 
reversed.  The  transducer  (11)  responds  cyclically  with  gradually  increasing 
amplitude  but  does  not  seem  to  have  a  permanent  offset.  The  transducer  at 
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(12)  shows  a  snail  cyclic  rasponsa  and  a  gradual  build  up  of  pora  pressure 
to  about  16  kPa  and  a  slow  decay. 

6.3  DETERMINATION  OF  THE  HELD.  ACCELERATION 

The  yield  accelerations  of  the  models  were  determined  from  results  of 
a  slope  stability  analysis  program,  the  program  was  based  on  the  method 
of  slices  outlined  by  Morgenstem  and  Price  (1965)  and  has  been  modified  to 
accommodate  constant  inclined  body  forces  due  to  lateral  acceleration 
during  an  earthquake.  First,  the  shear  strength  was  back-calculated  by 
coaqparing  the  results  of  a  computer  run  without  earthquake  loading  with 
the  results  of  the  static  failure  of  model  A.  This  model  was  consolidated 
at  80  g  after  a  series  of  four  earthquakes  and  the  speed  of  the  centrifuge 
was  gradually  increased  (over  about  four  minutes)  until  it  failed  at  122  g 
due  to  the  increased  centrifuge  acceleration  body  forces.  Hence  it  was 
known  that  the  factor  of  safety  was  1.0  at  122  g.  The  factor  of  safety 
of  the  embankment  at  80  g  was  then  122/80  *  1.53.  Frcm  this  known  factor 
of  safety  and  the  results  of  the  computer  analysis,  the  undrained  shear 
strength  could  be  calculated  to  be  19  kPa.  Results  of  repeated  computer 
runs  with  horizontal  acceleration  coefficients  of  0.1,  0.15,  0.2,  and  0.25 
and  the  known  shear  strength  are  plotted  in  figure  6.4  as  factor  of  safety 
versus  earthquake  acceleration  coefficient.  Frao  this  figure  it  was 
possible  to  determine  that  the  yield  acceleration  coefficient  (corresponding 
to  a  factor  of  safety  of  1.0)  was  ky  ■  0.2  . 

The  height  of  a  prototype  slope  is  proportional  to  the  centrifuge 
acceleration.  The  factor  of  safety  against  failure  of  an  undrained  cohesive 
slope  on  board  the  centrifuge  is  therefore  inversely  proportional  to  the 
centrifuge  acceleration  since  F.s.  Cu/yh.  So,  for  the  centrifuge  tests 
at  25%  greater  centrifuge  acceleration  (tests  AIV  and  CV)  the  factor  of 
safety  against  "static”  slip  is  1.53/1.25  ■  1.22.  The  factor  of  safety 
corresponding  to  any  lateral  earthquake  acceleration  coefficient  is  also 


reduced  by  25%.  Therefore,  the  lateral  acceleration  coefficient  that 
produced  a  factor  of  safety  of  1.25  for  the  6.9  m  slopes  will  produce  a  factor 
of  safety  of  1.0  for  the  8.6  m  high  slopes.  Osing  the  above  reasoning,  a 
relationship  between  factor  of  safety  and  lateral  acceleration  coefficient 
for  the  8.6  m  high  embankments  was  determined  by  scaling  down  the  relation 
for  6.9  m  high  slopes.  The  extrapolated  relationship  is  dotted  on  figure  6.4, 
the  yield  acceleration  can  be  picked  off  at  about  9%  of  the  centrifuge 
acceleration. 

6.4  PREDICTIONS  OF  DISPLACEMENTS  USING  THE  SLIDING  BLOCK  MODEL 

The  base  acceleration  data  from  FM  magnetic  tape  for  four  earthquakes 
was  digitized  and  stored  on  floppy  disc.  A  Fortran  program  was  written  to 
calculate  sliding  block  displacements  for  the  digitized  earthquakes.  The 
output  of  the  program  gives  relative  velocities  and  relative  displacements 
of  the  sliding  block.  The  results  of  these  calculations  are  summarized  in 
Table  III.  Plots  of  predicted  relative  velocity  and  displacement  with  time 
are  shown  in  figures  6.5  and  6.6.  These  plots  were  obtained  using  a  reduced 
shear  strength,  reasons  for  reducing  the  shear  strength  for  these  calculations 
will  be  explained  later.  These  plots  are  included  to  demonstrate  the  kind 
of  output  that  can  be  obtained  from  the  program. 

The  results  of  the  calculations  are  compared  with  measured  relative 
displacements  of  the  models  in  figures  6.5,  6.6.  and  6.7  and  in  Table  III. 

In  each  of  the  model  tests  there  was  either  one  or  two  displacement  transducers 
measuring  the  horizontal  displacement  of  the  slope  faces  at  mid  height.  The 
measured  displacements  that  are  compared  with  predictions  in  Table  III  are 
not  corrected  to  account  for  the  fact  that  horizontal,  near  surface 
displacements  are  not  identical  to  the  displacements  of  the  centre  of  gravity 
of  the  sliding  soil  mass.  Because  the  transducers  were  not  aligned 
exactly  in  the  direction  of  relative  motion,  and  were  not  anchored  at  the 
centre  of  gravity  of  the  sliding  soil,  they  may  have  recorded  displacements  of 


only  80%  to  90%  Of  the  relevant  displacement.  Because  this  correction  is 
indefinite ,  in  addition  to  being  small  compared  with  the  errors  in 
predictions,  it  was  preferable  to  present  the  unmodified  data  in  the  tables 
and  figures  of  this  report.  Figure  6.7  is  a  graphic  presentation  of  the  data  in 
table  III.  The  model  displacements  have  been  converted  to  prototype  displacements 
and  plotted  against  the  error  ratio:  predicted  displacement/measured  displacement. 

For  tests  BI,  CIII  and  DI,  the  predicted  displacements  were  all  less 
than  10%  of  the  measured  displacements.  Only  for  test  CV,  which  Involved 
much  larger  displacements ,  were  the  predictions  within  a  factor  of  two 
of  the  measured  displacements.  There  are  several  possible  causes  for  this 
large  discrepancy: 

(1)  Plastic  deformation  occurs  at  stresses  below  the  shear  strength. 

(2)  Dynamic  magnification  of  accelerations  In  the  embankment  produces 
larger  accelerations  and  Inertia  forces. 

(3)  •  The  soil  strain  softens  during  the  earthquake. 

(4)  The  sliding  block  model  restricts  shear  to  a  single  slip  surface 
whilst  in  the  centrifuge  models,  the  deformation  has  a  greater 
degree  of  freedom. 

The  fourth  of  the  listed  causes  requires  more  study  to  evaluate  its 
importance .  There  is  evidence  to  suggest  that  the  first  three  causes  did 
contribute  to  the  discrepancy  between  predictions  and  measurements.  These 
causes  will  be  discussed  below. 

(1)  Soil  is  an  elasto-plastic  work  hardening/softening  material.  It  does 
exhibit  plastic  strains  at  stresses  below  the  shear  strength.  Undrained 
shear  of  Kaolin  requires  2%  to  15%  shear  strain  to  mobilise  the  undrained 
shear  strength  in  monotonic  loading.  A  1%  shear  strain  might  mobilise  only 
about  50%  of  the  shear  strength  of  lightly  over-consolidated  Kaolin 
(Sketchley,  1973) .  A  1%  shear  strain  distributed  over  a  20  m  thick  shear 
band  could  produce  0.2  mm  of  displacement  per  pulse  of  acceleration,  or  2.0  mm 
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over  the  10  pulse  earthquake.  This  Is  of  the  same  order  as  the  displacements 
recorded  In  Tables  II  and  III.  An  Important  improvement  of  the  sliding 
block  model  would  be  to  account  for  the  non-linear  nature  of  the  soil  instead 
of  assuming  rigid  plastic  behaviour. 

One  method  of  accounting  for  yield  at  stresses  below  the  shear  strength 

would  be  to  use  a  non-linear  force  displacement  relation  in  the  sliding 

block  analysis.  A  simpler  method  would  be  to  use  a  reduced  but  constant 

strength  in  the  analysis.  This  involves  the  selection  of  an  arbitrary 

reduced  strength.  Some  calculations  were  conducted  using  a  reduced  strength 

of  90%  of  the  actual  strength.  The  results  of  these  calculations  are 

also  listed  in  Table  III,  and  are  used  for  the  plots  of  figure  6.5  and  6.6. 
Using  the  reduced  strength  improves  displacement  predictions  considerably, 
especially  for  the  earthquakes  that  produced  small  deformations.  However, 
for  test  CV,  the  reduced  strength  results  in  an  overestimation  of  displacements. 
These  results  suggest  that  the  effective  reduced  strength  depends  on  the 
magnitude  of  the  mobilised  stresses  and  displacements. 

(2)  Dynamic  magnification  of  accelerations  by  factors  between  1.25  and  2.0 
was  measured  in  every  earthquake  test.  Because  the  base  and  crest  motions 
often  had  different  frequency  content,  it  was  difficult  to  compare  their 
amplitudes  by  simply  picking  values  off  plots  such  as  figures  6.2  and  6.3. 

The  values  of  base  and  crest  acceleration  recorded  in  Table  II  were  taken 
from  filtered  acceleration  records  (filtered  by  eye  only)  -  the  high 
frequency  spikes  like  those  in  the  base  acceleration  traces  of  figure  6.2 
were  disregarded. 

Further,  more  objective  analysis  will  be  made  using  Fourier  Spectra 
and  comparisons  of  base  and  crest  velocities  rather  them  accelerations 
in  order  to  evaluate  the  importance  of  magnification. 

(3)  During  the  earthquakes,  pore  pressures  are  generated  and  the  soil  may 
be  partially  remoulded.  This  may  cause  a  reduction  in  the  yield  acceleration 
and  an  increased  rate  of  deformation  as  the  earthquake  progresses.  The 


-  19  - 


comparison  batman  predicted  and  measured  displacements  during  earthquake  DI 
in  figure  6.6  (using  a  reduced  strength)  demonstrates  a  phenomenon  that 
could  be  explained  by  strain  softening.  As  the  earthquake  progressed/  the 
amplitude  of  the  base  acceleration  actually  decreased.  In  the  displacement 
prediction,  this  was  manifested  as  a  reduction  in  the  rate  of  displacement 
during  the  latter  half  of  the  earthquake.  But, the  rate  of  measured 
displacements  actually  increased  as  the  earthquake  progressed. 

6.S  DELAYED  FAILURE  AFTER  EARTHQUAKE  CV 

The  most  catastrophic  of  the  observed  failures  caused  by  base  shaking 
was  that  after  test  CV.  This  model  sustained  large  deformations  (about  8  mm) 
during  the  earthquake.  Then,  about  3.5  seconds  after  the  earthquake  a 
collapse  occurred  during  which  about  70  mm  of  displacement  occurred.  The 
results  of  the  delayed  failure  are  shown  in  the  radiograph  of  figure  5.5c 
and  the  photographs  of  figures  5.7d  and  5.9.  Figure  6.8  shows  some  transducer 
readings  between  the  time  of  the  earthquake  and  the  "static"  collapse.  The 
top  trace  shows  the  base  acceleration  -  the  duration  of  the  earthquake  of 
about  0.17  s  occupies  about  5  mm  on  the  time  axis  near  the  beginning  of  the 
trace.  The  pore  pressure  responses  shown  in  the  next  two  traces  were  taken 
from  transducers  directly  beneath  the  crest.  The  first  was  at  an  initial 
depth  of  125  mm  from  the  crest  and  was  probably  in  the  sliding  mass  of  soil, 
but  very  near  to  the  shear  zone.  The  second  pore  pressure  was  measured  at  a 
depth  of  about  155  mm  from  the  crest  and  could  be  accurately  located  relative 
to  the  slip  plane  using  post  test  radiographs  (see  figure  5.5c.).  This 
second  transducer  was  approximately  20  mm  beneath  the  shear  zone  in  the 
part  of  the  embankment  that  did  not  "statically”  slip.  The  first  trace  shows 
a  small  negative  pore  pressure  generated  during  the  earthquake  gradually 
becoming  positive  as  the  collapse  occurred.  The  second  transducers  shows 

a  very  different  response  -  a  very  small  pore  pressure  generated  during  the 

■  t 

quake  and  then  a  gradual  build  up  of  excess  pore  pressure  to  about  9.0  kPa 
when  the  failure  occurred. 


Tha  horizontal  displacement  at  mid  height  and  the  crest  settlement 
traces  remain  steady  until  the  earthquake*  when  they  undergo  several 
millimetres  of  movement.  After  the  quake,  a  steady  creep  is  apparent  on 
both  traces  until  the  collapse  occurs  after  about  3.5  s. 

The  cause  effect  relationships  for  this  failure  cannot  be  understood 
until  further  study  cam  be  given  to  pore  pressure  data  in  the  after 
this  and  other  earthquakes.  However,  it  seems  likely  that  this  delayed 
failure  was  a  result  of  strain  softening  during  the  earthquake  deformation 
and/or  pore  water  pressures  generated  during  the  earthquake. 
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7.  CONCLUDING  REMARKS 

As  five  models  tested  gave  consistent,  repeatable  results.  Furthermore, 
models  of  different  scale  that  represent  the  same  prototype  agree  with  each 
other,  adding  support  to  the  centrifuge  scaling  laws.  The  Newmark  Sliding 
Block  Calculation  using  a  rigid  plastic  soil  strength  on  a  distinct  slip 
surface  consistently  under  predicts  measured  displacements.  Shear  zones  In 
clay  embankments  during  earthquakes  do  not  tend  to  form  distinct  rupture 
planes.  In  these  models  shear  zones  of  10-30  mm  thickness  were  observed.  Bence 
from  small  cyclic  displacements  of  less  than  1  mm  per  cycle  we  would  obtain 
less  than  3%  to  10%  strains.  These  small  strains  are  not  large  enough  to 
mobilise  the  full  undrained  shear,  strength,  but  can  result  in  permanent 
deformation.  Permanent  displacements  due  to  stresses  below  the  shear  strength 
must  be  accounted  for  in  order  to  improve  predictions.  Other  factors  which 
should  be  accounted  for  are  the  magnification  of  base  accelerations  due 
to  resonance  of  the  embankment,  and  strain  softening  during  the  course  of 
the  earthquake . 

The  collapse  of  model  CV  3.5  s  after  the  earthquake  resulted  in 
deformations  an  order  of  magnitude  larger  than  any  measured  during  the 
earthquake.  This  delayed  "static"  collapse  caused  by  strain  softening 
and/or  pore  pressure  generation  during  the  earthquake  may  be  the  most 
Important  consideration  when  designing  dams  to  resist  earthquakes. 

If  the  delayed  failure  is  due  to  pore  pressure  dissipation  or  redistribu¬ 
tion,  the  3.5  s  scales  to  about  two  hours  in  prototype  time. 
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TABLE  XI  Summary  of  Data 


Test  Acceleration 

(g) 


(Earthquake  Accn) 
(Centrifuge  Accn) 
Base  Crest 

(«)  («) 


Crest 

Settlement 

(mm) 


Horizontal  Displacements 


+Z  Side** 
(mm) 


-Z  side  ** 


AX 

80 

3.8 

7 

0.003 

0.003 

0.002 

All 

80 

13.5 

23.4 

0.218 

0.209 

0.166 

AXIX 

80 

24.5 

33.5 

2.025 

2.030 

1.58 

AXV 

100 

4.2 

8.4 

0.098 

0.055 

0.056 

BX 

80 

25.1 

31.2 

2.681 

1.985 

2.155 

BXX 

80 

2.1 

3.3 

0.001 

0.004 

0.003 

BXXI 

80 

7.9 

14.0 

0.009 

0.004 

0.009 

BIV 

80 

20.6 

31.2 

1.215 

0.835 

1.123 

3V 

80 

23.3 

34.0 

1.633 

1.240 

1.531 

Cl 

35.6 

2.5 

3.3 

0.012 

— 

0.012 

CII 

35.6 

9.5 

12.7 

0.023 

- 

0.00 

CIII 

35.6 

14.1 

18.6 

0.188 

- 

0.187 

CIV 

35.6 

17.9 

24.3 

1.224 

- 

1.22 

cv 

44.5 

17.7 

24.0 

9.0 

- 

7.6 

DX 

35.6 

21.8 

28.0 

3.254 

3.221 

2.205 

DXX 

35.6 

3.3 

7.0 

0.009 

0.0003 

0.004 

DUX 

35.6 

21.7 

32.1 

3.027 

4.71 

2.42 

El 

35.6 

.19 

24 

0.1 

0.059 

0.055 

EXX 

35.6 

21 

27 

0.45 

0.363 

0.257 

EXXI 

35.6 

23 

34 

0.69 

0.593 

0.347 

EIV 

35.6 

23 

34 

0.66 

0.601 

0.343 

EV 

35.6 

23 

34 

0.61 

0.554 

0.312 

**  +Z 

side  and  -Z 

side  refer  to 

faces  of 

the  model. 

The  +Z  side  Is 

the  face 

nearest  the  leading  edge  of  the  centrifuge. 


TABLE  III 


Summary  of  Sliding  Block  Predicted  Displacements  and  a  Comparison 
with  Measured  Horizontal  Displacements  at  mid-height  of  the  models. 


TEST 

Measured 

Displacements 

(mm) 

+Z  -Z 

Predicted  Displacements 
(mobilised  strength  *  Cu) 

ky*  +Z  -Z 

Displace-  :  Displace¬ 
ment  ment 

Predicted  Displacements 
(mobilised  strength  >0.9  Cu) 

ly*  +Z  -Z 

Displace-  Displace¬ 
ment  ment 

BI 

1.985 

2.155 

0.2 

0.162 

0.256 

0.15 

0.623 

0.837 

CIII 

0.187 

0.2 

0.001 

0.000 

0.15 

0.048 

0.085 

CV  • 

7.6 

0.09 

6.54 

5.49 

0.03 

23.5 

18.1 

DI 

3.221 

2.205 

0.2 

0.066 

0.208 

0.15 

0.622 

1.493 

»  Yield  Acceleration/Centrifuge  Acceleration 


\ 

I 


Emm 


IJil 


Schematic  view  of  Bumpy  Road  Apparatus 
looking  downward  into  centrifuge  pit.  The 
specimen  container  is  shown  in  its  swung  tip 
and  seated  position.  . 


Jtf/ml 


3.1b  Schematic  view  of  Bumpy  Road  Apparatus 

looking  radially  outward  from  centrifuge 
axis.  There  is  a  cutaway  view  through 
the  swinging  specimen  container. 


i 

v 
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Fig.  3.2b  A  tangential  view  over  the  top  of  the  centrifuge 
arm.  The  steel  track  is  held  on  the  wall  by 
brackets  that  are  enclosed  in  streamlining.  The 
nylon  wheel  is  being  held  clear  of  the  track  by  the 
pneumatic  jack. 


to  infinity 


Fig.  4.1  Relationship  between  normalised  displacement  and  the 
ratio  of  yield  and  base  accelerations.  (From  Kuttar, 
1979) . 


Fig.  4.3a 


ASSUMED  CONDITIONS  FOR  SENSITIVITY  ANALYSIS 

*  SLIDING  BLOCK  ANALYSIS 

*  2:1 ,  COHESIVE  SLOPE 

*  YIELD  ACCELERATION,  K^O.lg 

*  1  Hi  SINE  PULSE  . 

*  HORIZONTAL  BASE  VIBRATION 


K  g  sin(2n"t) 
nr 

Fig..  4.3b. 

PREDICTED  AND  ACTUAL  DISPLACEMENTS  IF  PREDICTIONS 

ARE- BASED.  QN.  A  5%  .QVEBESTIMATIQN..QF  SOIL -STRENGTH 


CASE 

Km 

Vk,, 

PREDICTED 

ACTUAL 

PERCENTAGE 

DISPLACEMENT 

DISPLACEMENT 

ERROR 

I 

0.11 

0.9 

1.0  mm 

17  mm 

1600.% 

II 

0.5 

0.2 

400  mm 

460  mm 

IB  % 

Fig.  5.3 

MODEL  B 

Dimensions  and  Instrumentation 
Vibration 


o  Pore  Pressure*  Transducers 
Accelerometers 

I - 1  r  LVDT  Displacement 

1 —  »  Transducers 


Lead  Threads 


i 


ELEVATION  VIEW 


10. 


Fig.  5.5d  Model  D  Earthquake  induced  distortion  apparent  near  the  base  of  the  specimen. 


5.6a  Model  B  photographed  in  flight  before  earthquakes  (SOg) 


* 


Fig  5.6b  Model  B  photographed  in  flight  after  all  earthquakes  (80g) 
Bolding  a  straight  edge  against  the  surface  markers  will 
reveal  distortion  of  the  embankment.  More  accurate 
measurements  on  the  film  measuring  machine  would  enable 
determination  of  volumetric  and  shear  strains  of  this 
section. 


Pig.  S.7a  Model  C  photographed  in  flight  at  35.6gf  before  any 
earthquake. 


Fig.  5.8b  Modal  D  photographed  in  flight  at  35.6g(  after 
all  3  earthquakes . 


Fig.  5.9b  Modal  C  after  removing  from  swinging  container 


PROTOTYPE  CREST  SETTLEMENT  (m 


8-6  m  high  embankments 


Pig.  6.1  Experimental  relation  between  prototype  crest 

settlement  and  the  maximum  earthquake  acceleration 
coefficient 
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Sressurtfi 


6.9m. embankment 
8.6m.  embankment 


ACCELERATION 
CENTRIFUGE  ACCELERATION 


Results  of  method  of  slices  slip  circle  analysis  on 
6.9  m  high  prototypes  and  extrapolated  curve  for  8.6  m 
high  prototypes. 


Figure  €.5  Input  base  acceleration,  calculated  relative  velocities  and 
relative  displacements  and  comparison  with  measured  displacements  for 
test  BZ.  For  calculations  the  strength  was  assumed  to  be  0.9  C  and 
the  corresponding  yield  acceleration  is  0.15(80  g)*»12  g. 
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Relative  Displacement 
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,  D  . .  .  .  Predicted  Displacement 

s  Measured  Displacement 


Relation  between  the  error  ratio  for  displacement 
predictions  and  total  prototype  displacement  for 
some  of  the  tests. 


Fig.  6.7* 


